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The dynamic phase diagram of vortex lattices driven in disorder is calculated in two and three 
dimensions. A modified Lindemann criterion for the fluctuations of the distance of neighboring 
vortices is used, which unifies previous analytic approaches to the equilibrium and non-equilibrium 
phase transitions. The temperature shifts of the dynamic melting and decoupling transitions are 
found to scale inversely proportional to large driving currents. A comparison with two-dimensional 
simulations shows that this phenomenological approach can provide quantitative estimate for the 
location of these transitions. 
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I. INTRODUCTION 



The prediction of nonequilibrium phase transitions in 
driven vortex lattices!] has triggered an extensive the- 
oretical study of periodic structures moving .through 
a random— epjironment, using both analyticalHu and 
numericalijiij appcaaches. The striking experiments on 
vortex transporttjEil that motivated and supported the- 
oretical efforts, provided convincing evidence for a gen- 
uine nonequilibrium phase transition between different 
driven states of the vortex lattice. InvestigationSpaf— the. 
I — V curves in Nb2Se samples and MoGe filmgiaM'ta 
revealed regimes of plastic and elastic flow of the vortex 
lattice separated by a.-characteristic current. Neutron 
scattering experimentscH on the driven vortex state in 
Nb2Se demonstrated a sharp increase in the density of 
structural defects of the vortex lattice in a certain cur- 
rent interval just above the depinning transition. At a 
larger characteristic current the defect density dropped 
significantly, suggesting a reordering of the vortex sys- 
tem. 

To understand the above experiments it is crucial to 
determine how quenched disorder affects the structure of 
the driven vortex system. The central idea of Ref. [I] was 
the suggestion that disorder, being frozen in the labora- 
tory frame, appears as a temporarily fluctuating force in 
the frame of the moving vortices and leads to an increased 
effective temperature of the vortex system. This effective 
temperature strongly depends on the drift velocity. Thus 
an equilibrium phase transition of the pure system (like 
melting) should have a counterpart in the system driven 
in a disordered environment and could be triggered by 
changing the drive at constant true temperature. For 
increasing drift velocities, disorder should be more and 
more washed out and the temperature of the dynamic 
transition should approach the transition temperature of 
the pure system from below. 

One related question of principal interest is to what 
extent disorder-induced features of the equilibrium phase 



diagram could be eliminated by applying a driving cur- 
rent. The most prominent example for such a feature 
is the first-order melting transition of the vortex lat- 
tice, which is experimentally found to end at a critical 
point H3 For decreasing strength of disorder the location 
of the transition moves to higher temperatures and the 
end point wanders to larger fields .c2l Therefore it seems 
plausible that with increasing sample purity this critical 
point continues to be shifted to larger fields and the theo- 
retically expected phase diagram, which does not display 
a critical point, is recovered. Since it is practically im- 
possible to obtain pinning-free samples, one might hope 
to observe the phase diagram of the ideal pure system 
by increasing a driving force in a given impure sample. 
The discussion of this scenario is one of the goals of the 
present paper. 

Already in the pure case one has to face the problem of 
how to capture the melting transition theoretically. Here 
a phenomenological approach basecL-Qn— the Lindemann 
criterion has been used successfullyOcS This criterion 
states that the static lattice melts as the mean squared 
thermal displacement of a vortex line becomes equal to a 
certain fraction of the lattice spacing, (u 2 (r, t)) ~ cf^a 2 ,. 
The number cl is called the Lindemann number and is 
usually of order unity, ao is the vortex spacing in the 
direction perpendicular to the magnetic field. In this 
conventional form the Lindemann criterion would suggest 
that in two dimensions (where (u 2 (r, t)) = oo for all finite 
temperatures) crystals would always be unstable to ther- 
mal fluctuations. Indeed, the long-range translational 
order is lost, but a quasi-long-range translational order 
and the topological order persist at low temperatures and 
vanish only above a finite melting temperature. Since dis- 
order in general reduces long-range order to quasLlong- 
range order between two an four dimensions it is 
necessary to modify the phenomenological criterion for 
a detection of the more subtle loss of topological order. 
This is achieved by using a slightly modified criterion for 
the relative displacement of two neighbored vortices, 
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W (b) = ([u(r + M)-u(r,i)] 5 



(1) 



Here b is a basis vector of the perfect lattice that sep- 
arates the undisplaced vortices. If one naively thinks 
of the lattice as being built by vortices connected by 
springs, then w is a measure for the typical stress of a 
bond ("spring") connecting neighboring vortices. Such 
bonds can be expected to "break" for 
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In a layered superconductor, where the vortex lines are 
actually composed of point-like vortices, not only, the 
melting transition but also a decoupling transitioned can 
exist. The evaluation of bond widths w for different ori- 
entations (in-plane of out-of-plane) provides additional 
insight into the anisotropy of vortex fluctuations and the 
location of these transitions. 

The extension ofr-fche Lindemann criterion on disorder- 
dominated systems^ enabled the description of disorder- 
induced static transitions between elastic and plastic 
glasses. In view of the loss of topological order the so- 
called "entanglement" transition from the Bragg glass to 
the vortex glass is equivalent to the melting of the pure 
system. Using the very same Lindemann criterion (||) 
with bond widths averaged over thermal and disorder- 
induced fluctuations, the location of this, entanglement 
transition has been determined recently.tETES The valid- 
ity of this criterion has even been derived within a self- 
consistent variational approachBj 

For driven vortex systems disorder-induced displace- 
ments were calculated in Ref. |l| within a naive large ve- 
locity expansion approach for the two-dimensional vortex 
lattice, and the nonequilibrium melting line was found 
for the case of strong disorder. Based on scaling ar- 
guments Balents and Fisher extended the concept of 
nonequilibrium, phase transitions over to charge density 
wave systems.H Giamarchi and Le DoussaB focused on 
the structure of the high velocity driven vortex phase and 
noticed that it retains some glassy features ofj-A^^^tic 
counterpart. The investigations that followedu&LJLl re- 
vealed further fundamental features of the driven phase, 
in particular, the fact that in a coarse-grained description 
(i.e., on a large enough spatial scales) the main effect of 
disorder on the driven periodic structure can be described 
as a random force. 

Most of these recent studies have been focused on the 
large-scale properties of the system in the elastic approxi- 
mation. Based on the properties of the topologically con- 
strained system it is in principle possible to include topo- 
logical defects like dislocations and to study the topologi- 
cal stability of the lattice systematically. In practice this 
approach is difficult to realize, in particular when the 
transition appears as collective effect of the dislocations, 
as in the example of three-dimensional melting. Encour- 
aged by the success of the Lindemann criterion in captur- 
ing the location of the phase transition in the pure case 
as well as in the static disordered case, we extend this 



approach in this paper to the location of nonequilibrium 
transitions of the vortex system. 

The organization of this paper is as follows: in Sec. || 
we outline the perturbative approach for vortex lattices 
driven through disorder and specify the generalized form 
of Lindemann criterion. In Sec. [II we evaluate the Lin- 
demann criterion for two-dimensional systems and per- 
form a quantitative comparison of the resulting phase di- 
agram with numerical simulations. In Sec. [V we derive 
consequences for the dynamic phase diagram of three- 
dimensional systems with a dynamic melting and decou- 
pling transition in 



II. PERTURBATIVE APPROACH 

We consider a vortex lattice in a steady driven state 
with velocity v. Vortices are labeled by their ideal posi- 
tion r in the comoving frame, whereas their actual posi- 
tion in the laboratory frame is R(r, t) = r + vt + u(r, i). 
The dynamics of the system is governed by the conven- 
tional equation for the vortex displacement field u = 
u(r,t): 



?7U = cV 2 u + f pin (R) 
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(3) 



where c represents the elastic constants (specified below 
for the vortex lattice) and V the lattice gradient. The 
thermal noise C couples the vortices to a heat bath of 
temperature T. The pinning force density f pia (r, t) — 
— W[R(r, t)] is related to the pinning potential V, which 
is supposed to be Gaussian distributed with a second mo- 
ment 



V(k)V(k / ) = A(k)<5(k + k') 



A(k) = A e 



(4a) 
(4b) 



The correlations decay on the scale of the coherence 
length £. We define k]_ = k%. + ky as the vector compo- 
nent perpendicular to the magnetic field. Eq. (||) refers 
to point-like vortices in a single layer (D = 2) or in a 
layered superconductor [D = 3). 

In the absence of disorder the response of the vortex 
lattice is different for longitudinal (L) and transverse (T) 
modes. Since the vortex lattice is almost incompressible, 
the longitudinal modes do not significantly contribute to 
many physical properties. However, as we see later on, in 
a vortex lattice driven through disorder these modes are 
important and have to be retained. The response func- 
tion G is determined by the elastic constants for com- 
pression Cn, shear cqq and tilt C44: 



G^(q,o;)=^Gf(q,a;)P^(q) , 



G p (q, uj) = {-irjuj + c p q\ + c A4 ql] 1 
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where p = L,T stands for a polarizations with c p = 
cn,C66 respectively, and g 2 = g 2 + q^. Wave vectors 
q are restricted to the first Brillouin zone (BZ) of the 
ideal lattice in contrast to k. 

To calculate the fluctuations w of the distance of neigh- 
boring vortices we treat disorder on the lowest level of 
perturbation theory by approximating f(R) = f(r + v£ + 
u) w f(r + vt), which is justified as long as w(h) < £ 2 . 
This condition is satisfied in the whole range of interest 
w(b) < c\a\ if clclq < £■ This restricts the validity of 
our approach to high magnetic fields not too much below 
the upper critical field. 

In this approximation the pinning force can be con- 
sidered as "external" force that does not depend on 
the response of the vortex lattice. The combination 
f(r, i) = f pm (r + vt) + C{ r >t) of this pinning force and 
the thermal noise has Gaussian correlations 



(/ Q (q, cj)f (q.', lj')) = ^(q, u)S(q + q')«(w + J), (6a) 
* Q/3 (q, to) = i9S aP + fcafc/3A(k)J(w + v • k), 
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i9 = 2r?T. 



The averaging includes the randomness of disorder and 
that of thermal noise. We denote k = Q + q with a re- 
ciprocal lattice vector (RLV) Q (being perpendicular to 
the magnetic field) and a vector q within the BZ. 

The displacements in response to the total external 
force f have correlations 



?i ct (q,w)u (3 (q', lj') = C Q/3 (q, w)6(q + q')6(uj + lj') (7a) 
C a/3 (q,u;) = G a7 (q,o;)* T «(q,w)G l 35(-q, -uj). 

From Eq. (^) we can immediately calculate the bond 
fluctuations 
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The thermal and pinning contributions to the force cor- 
relator (^|) generate two corresponding contributions to 
the bond width, w = w th + w pin , which can be considered 
as approximately independent. The thermal contribution 
has been studied extensively in the past in order to de- 
termine the melting transition of the vortex lattice in the 
absence of impurities. Here we focus our attention on the 
contribution due to pinning, which is 



pi "(b)«E / ^(q-b) 2 P>(k)A(k) 
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where J q = J jw^a and we abbreviate 
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rT(k) ^(k ± A qi 



r,2(vk)2 + ( C66 qi+ C4 4g z 2 ) 2 ' (10a) 



(10b) 



Before we proceed to a detailed evaluation of Eq. (||) 
in the following sections, we cite the result of lowest order 
perturbation theory for the macroscopic friction force F fr 
arising from the collective summation of the microscopic 
pinning forces. In order to achieve a drift velocity v one 
has to apply a driving force F(v) = ?yv + F fr (v) withB 

F* r = /ifc Q fc /3 fc 7 A(k)G /37 (k,-vk) (11a) 

= E / k ^ ■ k)r p (k)A(k) , (lib) 

which coincides for D = 2 with the early result of Schmid 
and Hauger.tj 

The disordered vortex lattice is characterized by two 
elementary length scales: the disorder correlation length 
£ and the "vortex spacing" a = ^/<f> /i? {B is the mag- 
netic induction and <i>o the flux quantum). In the tri- 
angular lattice the actual distance between neighboring 
vortices in a direction perpendicular to the magnetic field 

is a = J2/V3a - We assume c^ao < £ < ao which is 
realistic for high temperature superconductors at large 
fields. For the evaluation of the main formulae (j|) and 
( |lT| ) we will retain the reciprocal lattice structure in Q, 
but we approximate the BZ as spherical cylinder with 
bounds q^ < q*^ 2 = ^ja\ and \q z \ < q* z = n/d with the 
layer spacing d. This approximation preserves the area 
of the BZ. The nonlocality of cqq will be neglected. We 
choose the x-axis as the direction of the velocity, v = ux, 
and suppose the vortex lattice to move along one main di- 
rection of the hexagonal lattice —which are the directions 
of minimum energy dissipationO 



III. EVALUATION IN D = 2 

The general expressions for bond widths w and friction 
force F are immediately specialized to two dimensions 
by setting q z = 0. They are evaluated here in order 
to compare the location of the dynamic melting transi- 
tion according to the Lindemann criterion with numerical 
simulations. 



A. Thermal bond width 

To start with we consider the contribution of thermal 
fluctuations to the bond width. By apCpmparison with 
the Kosterlitz-Thouless melting theorycJ we are able to 
fix the Lindemann number cl- 

Thermal fluctuations lead to a displacement correla- 
tion 
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which implies a bond width 
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Here the contribution of longitudinal displacements can 
be safely neglected since typically c\\ 3> cqq. In D = 2 
there is only one bond length b = a. The actual—melt- 
ing temperature within Kosterlitz-Thouless theoryCJ (ne- 
glecting renormalization effects) is given by 
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Inserting this temperature into Eq. (13) one can estimate 
the Lindcmann parameter 
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0.03 



(15) 



This value is in good agreement with typical values that 
empirically describe the melting transition in high tem- 
perature superconductors. We stick to this value in all 
quantitative considerations below. 

As soon as we have determined the pinning bond 
widths we can estimate the location of the transition 
where bonds b become unstable according to the Lin- 
demann criterion w th + w pm = c 2 L a^. The transition 
temperature will be determined by 



T c — Tkt — 
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(16) 



as a function of the bond orientation and of velocity. 



B. Pinning bond width 

The evaluation of the disorder contributions ([))) to the 
bond width is more subtle than the calculation of the 
friction force. In particular, since the velocity selects a 
particular direction, which is chosen parallel to the x-axis 
(called IE-bonds), we have to distinguish bonds pointing 
parallel to the velocity and bonds that enclose an angle of 
60 degrees with the x-axis (called y-bonds). Both types 
of bonds have a length |b| = a. Due to the invariance of 
Eq. (^|) under a reflection of components of k, one can 
parameterize the bond width by two coefficients, 



(17) 



These coefficients w% m and wf, ln can be calculated ana- 
lytically for small and large velocities: 
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(18a) 
(18b) 



Upper and lower expressions hold for v 3> cu/a^ri and 
v <C ^c 6 g/aQ77 respectively. 

In the limit of large velocities the contributions to lead- 
ing order in v and in small £/ao come only from RLV 
with Q x = 0, since then the denominator of the response 
function becomes small, -q 2 v 2 k 2 + c 2 q 4 = ij 2 v 2 q 2 + c 2 q 4 , 
and gives largest weight to small q x . In the case of w%. m 
the elastic interaction is negligible above a characteris- 
tic velocity v ~ Cp/aorj. Then longitudinal and trans- 
verse modes are equally important and give contributions 
~ v~ 2 . For lug 111 vortices therefore respond as if they 
were independent particles. In the case of w^ ln there is a 
qualitative difference between longitudinal and transverse 
modes. Transverse modes again give only contributions 
~ v~ 2 , whereas now longitudinal modes give dominating 
contributions of order ~ Thus at large velocities 

?/-bonds are subject to much stronger fluctuations than 
x-bonds. 

For small velocities the elastic interaction dominates 
the response of vortices and for cu > c§q longitudinal 
modes are in general negligible compared to the trans- 
verse modes. The small velocity regime is reached when 
in the denominator rj 2 v 2 k 2 + c 2 q of the response func- 
tions the first term is typically small compared to the 
second one, i.e. below a velocity v ~ ^c v ja\r\. To the 
leading order in small v and £/<Xo, RLV of all directions 
contribute equally the the bond widths, which diverge 
~ ln(l/w). The anisotropy does not vanish at small veloc- 
ities, the prefactor of the logarithmic divergence depends 
(in subleading order in £/<Zo, arising from contributions 
of RLV with Q x = 0) on the bond orientation. In con- 
trast to the high velocity regime, x-bonds have stronger 
fluctuations than y-bonds at small velocities. 

These results will be discussed in more detail later in 



section HIE 



C. Friction force 



In most experiments and simulations the driving force 
rather than the velocity is the parameter controlling the 
drift. Therefore we wish to express the bond widths as a 
function of the force and to evaluate the transport char- 
acteristics v(F) from the friction force (pi]). Since we 
assume that vortices drift along a basic lattice direction, 
the friction force is also parallel to velocity, F fr || v || x. 



Eq. ( |11| ) can be evaluated analytically in the limits of 
large and small velocities. In both cases the main contri- 
butions to the sum over RLV's come from Q with all pos- 
sible orientations. Therefore the characteristic velocity 
separating the large and small velocity regimes for trans- 
verse and longitudinal modes are given by v ~ £c66/oq?7 
and v ~ ^cii/ayi] respectively. We find 
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for w > ^cn/a 2 ,/? , 
for i; < ^caa/alrj . 



(19a) 
(19b) 



Only for large velocities w ^> ^cn/a^rj longitudinal 
and transverse modes contribute equally to the fric- 
tion force. At velocities v <C ^Cn/a 2 ,^ the transverse 
modes dominate since en 3> cqq. Note that the fric- 
tion force enters the large velocity regime already for 
v 3> ^cn/dQfy, whereas the bond widths reach their cor- 
responding regime only for v 3> cn/aorj. 



D. Numerical evaluation 



agreement with the perturbative result is surprising for 
F > F c ss 0.04, where vortices move coherently (F c cor- 
responds to ft in Ref . . The regime of forces F < F c is 
beyond the validity range of our elastic approach. 

Fig. ^ displays the resulting phase diagram, compar- 
ing the result of the Lindemann criterion (full line) to the 
simulation (dots). Above this transition line all bonds are 
stable, below the transition line the Lindemann criterion 
breaks the breaking bonds. In the displayed range of 
small velocities ir-bonds are the least stable. Agreement 
between perturbation theory and simulations is given up 
to a factor of the order of unity, which is still quite fa- 
vorable in view of the conceptual simplicity of the Linde- 
mann approach and its sensitivity to changes e.g. of the 
Lindemann number. It is worthwhile to state that no fit 
parameters have been used in our numerical calculation. 



In order to illustrate the analytic results and to demon- 
strate the capability of the Lindemann approach to pro- 
vide even a quantitative estimate for the location of dy- 
namic phase transitions, we perform a numerical evalua- 
tion of the bond width and the friction force from Eqs. 
@ and (0). The results are then compared with the 
simulation data in Ref. [l]. 

For this purpose we specify the parameters as follows. 
In Ref. |l| the vortex spacing a was used as length scale and 
2deo = 2d(<E>o/47rA) 2 as energy scale (A is the penetration 
depth, d is the layer thickness). The time scale is set by 
i] = 1. From the melting temperature Tkt ~ 0.007 of the 
pure system we find Cqq » 



according to Eq. (14). 
The vortex interaction chosen in Ref. [j] decays on the 
penetration length estimated by A w a. The compression 
constant is obtained from en ~ (167rA 2 /a 2 )c66 ~ 50c66-til 
We furthermore identify £ = r p = 0.2 and refer to 
the data sets with a pinning strength A — 0.006, a 
number of pinning centers N p = 10 4 , and a number 
of vortices N v — 400. Then the disorder strength is 
A = lu/al = {N p /N v al)[TTA{S,/a ) 2 ] 2 sa 1.42 • 10~ 5 . 

Fig. [l] shows the pinning bond widths for x-bonds (full 
line) and y-bonds (long dashed line). The short dashed 
lines are a guide to the eyes representing an asymptotic 



decay 



for x-bonds and 



for y-bonds. This 



asymptotic regime is reached only for v > 100. At small 
velocities v < 0.01 both bond widths diverge logarithmi- 
cally. For large velocities, y-bonds have stronger fluctu- 
ations than x-bonds, as opposed to small velocities. 

The friction force obtained from Eq. ( ^l|) is shown in 
Fig. H by the full line. The dashed line displays aw -1 
dependence, which is realized by F for v > 100. At 
small velocities, v < 0.01, F fr saturates at a finite value 
in accordance with Eq. (|19|d). The resulting transport 
characteristics v(F) is shown in Fig. The dashed line 
is the characteristics in the absence of pinning, which is 
shifted to the full line by the presence of disorder. This 
shift is practically constant in the small velocity regime. 
The dots represent the simulation data of Ref. [I] for the 
lowest temperature (T = 0.001) considered there. The 



E. Discussion 



At this point we rest to discuss some specific properties 



of our results (18a) and (18b) for the bond widths 



The first property is the strong anisotropy at large ve- 
locity, where the width of a;-bonds scales like ~ 1/v 2 
independent of the elastic constants, whereas the y-bond 
width scales like ~ 1/wcn. This is true only for largest 
velocities v cn/aor), whereas for v <C cn/ao^y the 
isotropy is essentially restored, i.e. y-bonds have qualita- 
tively the same velocity dependence as x-bonds but the 
prefactors are still different. At large v the physical origin 
of the anisotropy can be understood in a simple picture 
where vortices are considered as independent particles. 
Vortices neighbored in x-direction follow the same paths 
and are exposed to the same pinning forces. They expe- 
rience the same force, but with a delay time At = a/v. 
Even in the absence of vortex interactions such purely 
time shifted forces give rise to bond fluctuations of a fi- 
nite width only. Therefore the elastic interaction is ir- 
relevant and these bond fluctuations are independent of 
elastic constants and decay proportional to (A£) 2 tx v~ 2 . 
Vortices neighbored in directions not parallel to velocity 
always move on different trajectories and are exposed to 
essentially uncorrelated pinning forces. In the absence of 
interactions their typical relative distance would increase 
without limits in the direction perpendicular to v as in a 
diffusion process. In the lattice such a diffusive motion is 
prevented by the vortex interactions. Since their relative 
distance fluctuates mainly in y-direction, which is almost 
parallel to their distance, it is the compression modulus 
rather than the shear modulus that confines the bond 
fluctuations. 

The anisotropy at large v can be related to 
anisotropy of the Larkin domain in the driven lattice.1 
This domain is much longer in the direction parallel to 
the velocity than in the other directions, which means 
that relative displacements grow faster in direction per- 
pendicular to v than parallel to v. Our conclusion that 



5 



j/-bonds are less stable than x-bonds at large velocities 
is in agreement with a smectic structure of the drift- 
ing system,ElQ~i!l where vortices move in decoupled chains 
aligned parallel to v. 

As second distinct feature of the bond fluctuations, the 
divergence of the bond widths at small velocities deserves 
some explanation. This divergence, which is in general 
found for D < 2, implies that at zero velocity and for ar- 
bitrarily weak disorder all bonds are broken and the vor- 
tex lattice is destroyed, i.e. that the structure factor re- 
sembles that of a liquid. This result of the Lindemann cri- 
terion coincides surprisingly with that of more elaborate 
methods, e.g. renormalization group methods rWA which 
evaluate the large-scale correlations of the displacement 
field. The divergence of the bond width in our approach 
arises from the integration over small momenta q, i.e. 
from large scale fluctuations. From a principal point of 
view our perturbative treatment of disorder is not valid 
at largest scales, where higher order effects can no longer 
be neglected. Therefore the correctly found divergence 
of the bond width has to be considered as lucky circum- 
stance. 

Although the phenomenological Lindemann criterion 
is able to capture the position of a given transition, it 
can neither proof the existence of a transition nor can it 
tell us something about the nature of a transition. Al- 
though in 2D a phase transition can be observed in the 
density of lattice defects and, as a consequence- in spe- 
cific features of the transport characteristics Jl3£3 scaling 
argum£n±a3 and numerical calculations of the structure 
factonilJEj indicate that on large length scales free defects 
should exist, i.e. that the topological order of the lattice 
is lost. However, at present a rigorous characterization of 
the large-scale properties and examination of the effects 
of defects is still missing. Even if a true phase transi- 
tion (like "solid-to-fluid" ) related to a qualitative change 
of the large correlations might be absent in 2D case, we 
expect a transition (like "liquid-to-gas" ) or at least a pro- 
nounced crossover at the location of the dynamic phase 
"transition" located by the Lindemann criterion, which 
in its generalized form used here probes small scale corre- 
lations. While these displacement fluctuations on small 
scales are related to the rate of generation of (initially 
bound) dislocation pairs, large scale fluctuations of the 
displacement field are relevant to decide whether dislo- 
cations remain bound in pairs or dissociate into free dis- 
locations that destroy the topological order. 

It is instructive to draw a comparison between the 
bond widths, which we examine here, and the shaking 
temperature T s h introduced in Ref. |l|. T s h was defined 
from the correlator of the pinning force experienced by a 
single particle. It is therefore independent of elastic con- 
stants. If we compare the thermal and disorder contri- 
butions to the bond width, the latter could be expressed 
in terms of an effective "bond shaking temperature" 



This bond shaking temperature differs from T s h in sev- 
eral respects. First of all it depends on the orientation 
of the bonds under consideration. At large velocities T s h 



m com- 



and T for y-bonds have the dependence • 
mon. However they differ in the prefactor, which is inde- 
pendent of the elastic constants in the former case (since 
T s h characterizes a single particle) , but contains a prefac- 
tor ~ cgg/cn in the latter case (since T bsh characterizes 
a relative displacement response). T bsh is more similar to 
the coherent shaking temperature T^ h of Ref. [j], which 



was found to decay 



in agreement with T^h for 



x-bonds (and even y-bonds in an incompressible lattice). 

Speaking about effective temperatures, the attribute 
"shaking" should not be taken too literally, since even 
for v = the bond width w pm is finite and static, i.e. it 
has not temporal fluctuations. However, a snapshot of 
the displacements on small scales looks like that of a sys- 
tem without disorder but with a temperature increased 
by T hsh . Since an effective temperature cannot be defined 
uniquely in a nonequilibrium situation anyway, we keep 
in the following the notion of bond fluctuations. Since 
these fluctuations are related to the fluctuations of the 
Peach-Koehler force acting on dislocations, we believe 
that the bond fluctuations are an appropriate measure 
for the relevance of dislocation (at least on not too large 
scales, see discussion above). The order-of- magnitude 
agreement on the location of the transition between the 
Lindemann criterion and simulation supports this pic- 
ture. 



IV. EVALUATION IN 3D 

In a layered three-dimensional superconductor the 
physics is even more rich than in two dimensions: be- 
sides the melting transition, where the structural order 
of the vortex lattice gets lost and its elastic moduli get 
renormalized to zero, an additional decoupling transi- 
tion, where the conductivity perpendicular to the lay- 
ers becomes ohmic and the effective Josephson coupling 
between the layers gets lost, can occur. We evaluate 
the bond fluctuations for layered superconductors in the 
driven disordered case and discuss the implications for 
the nonequilibrium counterpart of the equilibrium melt- 
ing and decoupling transitions. 

In bulk superconductors, in particular for anisotropic 
high-temperature superconductors, the elastic properties 
of the vortex lattice are somewhat intricate due to the 
dispersion of the elastic constants (niamly of C44 and cu). 
For completeness we cite the valuestanH 
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C66 



1 +7 2 A 2 g 2 



al 1 + 7 2 A 2 g 2 + \2 q 2 i + AV 



T bsh (b) = ^l w(h) 
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We define q* ± = Vln/ao, q* z = n/d, and e = ($o/47rA) 2 . 
The tilt modulus is composed of a nonlocal contribution, 
to the tilt energy and a "single-vortex" contribution.EZl 
We are interested in the range of not too small fields, 
where A > clq and an exponential decay of the elastic 
constants with the vortex density can be neglected. 



A. Thermal bond widths 

It is straightforward to evaluate the bond fluctuations 
due to thermal fluctuations and in the absence of dis- 
order. Although purely thermal fluctuations have al- 
ready been Willi examined within the usual Lindemann 
approaches E3t3, we calculate the bond fluctuations in 
order to gain additional insight into the anisotropic na- 
ture of the fluctuations. Neglecting again the contribu- 
tions of compression modes, we find 



,.th _ 
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where 







2 _ C66<7j 



C44C 



B 

B rT 



is the magnetic field in units of the crossover fief 



7T$0C44 

4d 2 c 66 



(23) 



(24) 



For the transparency of our arguments and for qualita- 
tive purposes we have suppressed the explicit nonlocality 
of the elastic constants. The nonlocality can be restored 
by using effective values of the dispersive elastic constants 
( pi] ) evaluated at the length scales that give the largest 
contributions to the integrals. In D = 3 bond fluctu- 
ations are always dominated by small-scale fluctuations, 
since even in the presence of disorder and for small veloc- 
ities one never encounters a divergence arising from large 



scales (unlike in D < 2). Therefore we use in our fol- 
lowing estimates elastic constants evaluated for q± = q\ 
and \q z \ = q*. In general this approximation underes- 
timate the stiffness of the lattice but it is qualitatively 
valid as long as the dispersion (anisotropy of the elastic 
constants) is not too pronounced. 



Let us now locate the lines where w 



th 



r 2 n 2 
c l "0 



for all 



bond types a = x,y,z. These lines can be described 
by T a = T a (B) in the (T,B) plane. To be specific, we 
consider in the following moderately anisotropic systems 
with 7 d < A (like for YBCO, and even BSCCO lies just 
at the border) at not too small fields (do hS-A)- For such 
systems one finds at q± = q^ and q z = q* z ^\ 
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(25a) 
(25b) 

(26a) 
(26b) 



(22b) max(-u/ : 



The relative strength of fluctuations of in-plane- and z- 
bonds is different for small and large fields. At large fields 
w^/ciq is field-independent, whereas w* h / a o ~ 1/lnB 
increases without bounds, i.e. w* h <C u>' h . Thus for in- 
creasing fields T x saturates at a constant value, whereas 
T z vanishes. At small fields the situation is reversed: 

u4 h > u>* h because w* h / a o ~ B and w ^ h /«o ~ 51/2 ■ 

The bond widths provide more information than the 
typical displacement (u 2 ) 1 / 2 evaluated by the usual Lin- 
demann criterion. However, they are related by (u 2 ) w 

) for the 



th „,,th 



w. 



J ). In our Lindemann criterion 



different types of phase transitions we use the same Lin- 
demann number cr,. In principle these values might 
be different, since z-bonds and in-plane bonds are not 
not equivalent from the point of view of lattice symme- 
tries. But a use of different Lindemann numbers would 
lead only to a quantitative shift of the phase boundaries, 
which shall not be our main concern here. 

In order to associate the bond widths to possible melt- 
ing and decoupling transitions, it is worth recalling the 
main features of these transitions. Since for anisotropics 
7> 1 the interaction between vortices in the same layer 
is much weaker than the interaction of vortices in differ- 
ent layers one might in principle expect the possibility 
that melting occurs as a two-step process. In a first step 
the positional correlations of vortices in different layers 
could get lost while the correlations within the layers are 
preserved (at least as quasi-long-range order). This re- 
maining two-dimensional order within the layers could 
then become short-ranged at a second transition at a 
higher temperature. However, this scenario is probably 
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not realized, as one can conclude from an analogy to lay- 
ered XY-raodeh. On symmetry grounds their phase dia- 
gram should be topologically equivalent to that of layered 
crystals. It was showne3 that these models have only a 
single transition where order is destroyed in all directions. 
Although the relative fluctuations between the layers can 
become arbitrarily large for weak coupling between the 
layers, which is the case for layered superconductors at 
large fields, the transition temperature does not drop 
below that of the two-dimensional system. In the Lin- 
demann approach melting has therefore to be identified 
with the breaking of in-plane bonds. 

The breaking of z-bonds therefore describes a differ- 
ent transition, which is not related to the loss of crys- 
talline order. Two possible transitions have been con- 
sidered recently: the decoupling transition and the su- 
persolid transitioBp-T-hc decoupling transition of layered 
superconductorsEZlocJ describes the loss of phase coher- 
ence in the superconducting order parameter between dif- 
ferent layers and in general does not coincide with the 
melting transition. The location of this transition can be 
estimated by a Lindemann-type criterion for the phase 
difference between neighboring layers. Although the vor- 
tex displacements are the main source of such phase fluc- 
tuations, a Lindemann criterion for phase differences is 
in general not equivalent to a Lindemann criterion for 
vortex displacements, since a local vortex displacement 
leads to a nonlocal phase distortion. However, at fields 
below the crossover field B cr , both types of Lindemann 
criteria actually are equivalent and the decoupling, tran- 
sition is captured just by breaking of z-bonds.E3c3 At 
large fields the breaking of z-bonds describes the loca- 
tion of the transition from a solid to a supersolid vortex 
crystal due to a proliferation of_ vacancy and interstitial 
lines in the vortex line crystalo 



B. Pinning bond width 

After the brief review on the effect of thermal fluctua- 
tions on bond widths and the possible interpretations of 
bond-breaking we now focus on our main issue, the con- 
tribution of pinning to the bond widths and the question 
how the above equilibrium phase diagram of the pure 
evolves into a noncquilibrium phase diagram of the dis- 
ordered system. 

Only very recently the equilibrium phase diagram of 
the disordered system (v = 0) has received some atten- 
tion (see Refs. 32-38). In particular the breaking of in- 
plane bonds has been used as indication for the transi- 
tion from a topologically ordered phase (Bragg-glass) to 
a topologically disordered phase (vortex glass). For small 
magnetic fields the quantitative evaluation of the bond 
widths turned out to be quite a subtle issue because of 
the strong dispersion of the elastic constants in layered 
materials and the breakdown of perturbation theory in 
the low-field regime (£ <§C C£a ). Since we are eventu- 



ally interested in the dynamic effects of disorder on the 
dynamics, we do not attempt to reproduce the static re- 
sults. However, it is worthwhile to point out that in 
D = 3 the disorder induced bond fluctuations are finite 
for all bond types. This finitcness is found even at zero 
velocity and even if the bond widths are calculated from 
Eq. (^|) within naive perturbation theory that overesti- 
mates the effect of disorder. For weak enough disorder all 
bonds are stable and a topologically ordered phase can 
persist. For low temperatures the main effect of weak 
disorder is to reduce the long-range translational order 
to a quasi-long-range order. 

At sufficiently large velocities the perturbative ap- 
proach, Eq. (||), can always be used to calculate the shift 
AT a of the bond breaking temperatures T a with respect 
to the values T° in the absence of disorder. For large 
velocities these temperature shifts tend to zero and the 
use of the lowest order perturbation theory is justified at 
least as long as the disorder contributions to the bond 
widths satisfy w pin < £ 2 . 

The pinning bond widths show a complicated veloc- 
ity dependence with several regimes, separated by char- 
acteristic velocities related to the different elastic con- 
stants. We focus here exclusively on the high-velocity 
regime v 3> Cu/aof) for which we estimate the pinning 
contribution to the bond widths from Eq. (^): 



a 2 A 



w 



8Traodt; 3 r] 2 v 2 
q 2 A 
' ~ 4\/27r^ 3 ^wcii 
7ra A 



pin 



„P ln 



1 6 v 2-K^ A r}v^/ci i C44 



(27a) 
(27b) 
(27c) 



where we used again cng| 2 3> c^q* . 

From Eq. ( p7j ) one can immediately draw the impor- 
tant conclusion that y- and z-bonds, for which ui pin de- 
cays like are affected by disorder at large veloci- 



ties much stronger than x-bonds with w^." 



As 



in D = 2, the relative strength of y-bond fluctuations 
in comparison to a;-bond fluctuations is in agreement 
with the-.prediction of a smectic-like flow of the vortex 
lattice.H3 

In order to extract the net field dependence of the bond 
widths it is important to remember the implicit depen- 
dences of the parameters. In particular Ao ~ Oq 4 since A 
is the correlator of the pinning energy density in Eq. (Q). 
For the same reason 77 ~ <Jq 2 and also all elastic constants 



carry a dominant implicit field dependence 



Thus 



the relative fluctuations w^/oq ~ a are essentially field 

independent, whereas w% in /a§ ~ w £ m / a o ~ a Q 1 increase 
for large fields. 

To visualize the velocity dependence of the dynamic 
transitions it is instructive to translate the pinning bond 
widths (27) of the driven vortex lattice into shifts AT X of 
the bond-breaking temperatures rtgq) of the pure system. 
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From + u>P m = (? L a^ the bond breaking temperature 
is reduced by AT a = —T a w^ m /c 2 L al. These shifts are 



AT, 



AT,, 



AT, 



£qAo 



5 



-T 



-T, 



4:Trci£, 2 r]fv 2 \H c2 
eoA 



1/2 



7re A 



1/2 



(28a) 
(28b) 

(28c) 



To make the net field dependence of bond widths explicit, 
we have introduced a dimensionless disorder strength Ao 
by A = d(£eo/a§) 2 A and the friction coefficient per 
line iji = a^rj. Furthermore we have used again the val- 
ues of the elastic constants for q± = q*j_ and q z = q*, 
where en « 4c66 ~ £o/ a Oj and iJ C 2 = $o/27r£ 2 . In terms 
of these parameters the high- velocity regime is restricted 
by r\\v > e /a . 

In this high velocity regime y-bonds are less stable than 
x-bonds over the whole field range. Dynamic melting is 
thus always driven by the breaking of y-bonds. In the 
absence of disorder the z-bonds were more (less) sta- 
ble than in-plane bonds for low (high) fields compared 
to B cr . This relation is preserved in the driven case, 
since the bond breaking temperature of in-plane bonds 
is shifted more (less) than that of the z-bonds at low 
(high) fields. Thus, the topology of the phase diagram 
of the driven vortex lattice remains unchanged at large 
velocity, as one could expect naively. The high-velocity 
expansion does not provide evidence for a dynamic shift 
of the crossover field. 



V. DISCUSSION AND CONCLUSIONS 

To discuss the actual observability of the high- velocity 
regime it is important to compare the above velocity 
range to the depairing velocity of Cooper pairs. From 
the depairing current (see e.g. Ref. fll|) the depairing ve- 
locity can be estimated as 77; impair ~ eo/£- Thus even if £ 
is not much smaller than ao (as required by the validity 
of the perturbative approach) the high-velocity regime 
covers a substantial current range below the depairing 
current and should be experimentally accessible. As long 
as disorder is weak (A <C 1) the location of the dynamic 
melting and decoupling transitions is close to the loca- 
tion of the equilibrium transition of the pure system, i.e. 
AT Q « T a . 

In this high- velocity regime dynamic melting is always 
driven by the breaking of y-bonds. The dynamic shifts 
of both melting and decoupling transition temperatures 
scale like AT a ^ji -1 . This dependence is in a greem ent 
with experimentala and numerical simulationsBlli and 
the original concept of the shaking temperature. El We 
have extended this original approach, which was focused 
on the dynamic response of a single vortex and on the 



response of the lattice at large scales, on the fluctuations 
of bonds. Thereby we were able to put the character- 
ization of the noncquilibrium transitions on a common 
basis with the Lindemann approaches to the equilibrium 
counterparts. 

At smaller drive the velocity dependence of the bond 
widths and the temperature shifts will be in general 
weaker than at high drive. In two dimensions the velocity 
dependence becomes logarithmic at small velocities and 
the crossover between different regimes has been calcu- 
lated numerically (see Fig. |l|). In three dimensions there 
are even more dynamic regimes between the asymptotic 
regimes of zero and large velocity because of the complex 
elastic properties of the vortex lattice. 

However, from Eq. (28) we can conclude that for weak 
disorder (Ao <C 1) the scaling of the high- velocity regime 
should be observable in a considerable velocity range be- 
low depairing. In this case, where the temperature shifts 
are small (AT Q <C T a ), it should be possible to observe 
a dynamic shift of the phase boundaries experimentally. 
To the best of our knowledge, such dynamic shifts have 
not been reported. We naturally expect a dynamic shift 
to be most pronounced in the vicinity of the critical point 
of the melting line since the main effect of a driving cur- 
rent is to reduce the effective strength of disorder and this 
point has been found to be qmte sensitive to the strength 
of disorder in equilibrium .Ea In particular the location 
of the critical point will be shifted to larger fields by a 
driving current. 

The anisotropy of bond fluctuations found abo 
consistent with the picture of smectic vortex Ao 
due to the dominant proliferation of dislocations with 
Burger's vectors parallel to the drift velocity providing 
the actual mechanism for bond breaking. We expect the 
anisotropy of bond fluctuations to have immediate conse- 
quences for the structure of the coherently flowing vortex 
lattice. In particular even the Bravais basis of the vor- 
tex lattice will display this anisotropy. In the presence 
of disorder the vortex system will show an expansion 
of its lattice constants similar to a thermal expansion. 
From the anisotropy of the bond fluctuations one has to 
expect that this expansion will be larger in ^-direction 
than in x-direction. It would be interesting to analyze 
neutron scattering and Bitter decoration experiments on 
the driven vortex lattice in view of a dynamic shift of the 
location of the (quasi)-Bragg peaks. 
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FIG. 1. Plot of pinning bond widths w^, ln (full line) and 
Wy ln (long dashed line) calculated numerically from Eq. ([]) 
for the parameters specified in the text. The short dashed 
lines represent dependences w pln (w) ~ w _1 ,w -2 , which are 
realized at high v. 
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FIG. 2. Plot of friction force (full line) calculated numer- 
ically from Eq. ( pd| ) for the parameters specified in the text. 
The dashed line represent the dependence F lr (v) ~ u _1 at 
high v. 
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FIG. 3. Plot of the transport characteristic resulting from 
Eq. ([ll]), full line, and the simulation data of Ref. [l] (dashed 
line). 
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FIG. 4. Plot of critical force for the transition from co- 
herent to incoherent vortex motion. The full line is our result 
from the Lindemann criterion, the dots are simulation results 
from Ref. 
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